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ABSTRACT 

We present 3D hydrodynamical simulations of ram pressure stripping of massive disc 
galaxies in clusters. Studies of galaxies that move face-on have predicted that in such 
a geometry the galaxy can lose a substantial amount of its interstellar medium. But 
only a small fraction of galaxies is moving face-on. Therefore, in this work we focus 
on a systematic study of the effect of the inclination angle between the direction of 
motion and the galaxy’s rotation axis. 

In agreement with some previous works, we find that the inclination angle does not 
play a major role for the mass loss as long as the galaxy is not moving close to edge-on 
(inclination angle < 60°). We can predict this behaviour by extending Gunn & Gott’s 
estimate of the stripping radius, which is valid for face-on geometries, to moderate 
inclinations. 

The inclination plays a role as long as the ram pressure is comparable to pres¬ 
sures in the galactic plane, which can span two orders of magnitude. For very strong 
ram pressures, the disc will be stripped completely, and for very weak ram pressures, 
mass loss is negligible independent of inclination. We show that in non-edge-on geome¬ 
tries the stripping proceeds remarkably similar. A major difference between different 
inclinations is the degree of asymmetry introduced in the remaining gas disc. 

We demonstrate that the tail of gas stripped from the galaxy does not necessar¬ 
ily point in a direction opposite to the galaxy’s direction of motion. Therefore, the 
observation of a galaxy’s gas tail may be misleading about the galaxy’s direction of 
motion. 

Key words: galaxies: spiral - galaxies: evolution - galaxies: ISM - galaxies clusters: 
general 


1 INTRODUCTION 


Galaxies populate different environments, ranging from iso¬ 
lated field galaxies to dense galaxy clusters. Comparative 
observations between field and cluster galaxies reveal that, 
in addition to internal processes, the environment plays 
an important role in galaxy evolution. Disc galaxies are 
especially affected. Many well-known differences between 
cluster and field disc galaxy populations are attributed 
to a smaller gas content of cluster members. Such differ¬ 
ences are the higher fraction of red galaxies and early- 

S galaxies in clusters (e.g. iGoto et all 1200. 'll : IPimbbletl 
as well as the Hl-deficiency of cluster di sc galax- 
la.vatte et IDEMST olancs et al. 2Q01) and_J,he 

reduce d star formation ra tes j Koopmann fc Kennevlll998l 

2004^0 iKoonmann et, a.fll2005lh 
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Several processes have been proposed to explain these 
features. One idea are tidal interactions between cluster 
members. Although the interaction timescales in clusters are 
rather short due to the high velocity dispersion of cluster 
galaxies, Moore et al. (1996,1998,1999) have shown that re¬ 
peated short interactions (harassment) can have a s ubstan¬ 
tial influence on cluster galaxies (see also review bv iMiho^ 
120041 and references therein). Another process is ram pres¬ 
sure stripping, which is the interaction between the galactic 
gas disc and the intracluster medium (ICM). As galaxies 
move through the ICM, the ram pressure is expected to 
push out (parts of) their interstellar medium (ISM), pro¬ 
vided the ram pressure is strong enough. This process does 
not affect the dynamics of the stellar and dark matter com¬ 
ponents of the galaxy. If ram pressure stripping plays a role, 
one expects to find galaxies with undisturbed stellar discs 
and distorted and truncated gaseou s discs. Su c h examp les 
have been o bserved e.g. N GC 4522 Jkennev fe Koop mannl 
! 199)1, [2001; Kennev et al 1 l2004llVollmer et aDUOMl . NGC 
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4548 iVollmer et alJ IlflflflF) and NGC 4848 llVollmer et alJ 
1200 ill . 


The ICM-ISM interaction is a complex process in¬ 
fluenced by many parameters. D ifferent aspects hav e 
been studied b y sever al grou ps. ISchulz fe Struck l|20(nj), 
iMa.rcolini et a.l] d20n.‘4l - > and iRoedigei^^Henslen d2005tl 


showed that the ICM-ISM interaction is a multi-stage pro¬ 
cess. The most important phases are the instantaneous strip¬ 
ping, on a timescale of a few 10 Myr, an intermediate phase, 
on a timescale of up to a few 100 Myr, and a continuous strip¬ 
ping phase that, in principle, could continue until all gas is 
lost from the galaxy. Instantaneous stripping is the immedi¬ 
ate consequence of the ram pressure dislocating the gas disc. 
However, the gas does not become unbound from the galac¬ 
tic potential immediately. It takes the i ntermedia te phase 
until all displaced gas is truly lost. ISchulz fc Struck! ( 2001 1 
described this as a “hang-up” phase where gas lingers behind 
the disc. The continuous stripping is caused by the Kelvin- 
Helmholtz (KH) instability induced by the ICM wind flow¬ 
ing over the surface of the remaining gas disc, leading to a 
slow but continuous gas loss. This process is also termed tur¬ 
bulen t /viscous stripping (see e.g. iNnlsenlll 982l : lOnilis et a.lJ 

I2000I1 . 


A list of parameters that influence ram pressure strip¬ 
ping includes the ICM wind density, Picm, wind velocity, 
uicm, the inclination angle, i, between the disc’s rotation 
axis and the wind direction, the overall galaxy mass (its 
potential depth) and its gas mass, the variability o f th e 
ICM wind, the structure of the ISM gas disc. lCunn fc Gotti 
dl972T) proposed that for galaxies moving face-on through 
the ICM the success or failure of ram pressure stripping 
can be predicted by comparing the ram pressure with the 
galactic gravitational restoring force per unit area. If this 
comparison is done for each galactic radius, a stripping ra¬ 
dius can be estimated. In general, numerical simulations 
agree with this analytical estimate. For inclined geome- 
tries, no clear pictu re has yet emerged. In SPH simulations, 
lAbadi et al.l j l flflilti found a good agreement with the ana¬ 
lytical estimate for face-on cases. For the few inclined cases 
they studied, their model galaxies lost significantly less gas. 
However, this result may be biase d by their short runtime 
of ~ 100 Myr. Using an SPH code. ISchulz Kr. Strnc.kl d200lh 
studied the phases of the ICM-ISM interaction in detail. 
They focused on face-on geometries, but also performed a 
few simulations with inclined galaxies. They found that in¬ 
clined galaxies are stripped on longer timescales. They also 
studied the loss of angular momentum and found that in¬ 
clined galaxies are “annealed” by a stronger loss of angu¬ 
lar momentum. Hence, the galaxies co ntract an d thu s be - 
come more resistant against stripping. iQuilis et alJ ( 2000 ) 
used a hydrodynamical grid code to study, both, the influ¬ 
ence of the inclination angle and an inhomogeneous gas disc 
structure. In their work, they concentrated on rather high 
ram pressures, for which they found the inclination angle 
to play a minor role as long as th e galaxy is no t moving 
strictly edge-on. So far the work of lVollmer et al.l feOOlIl is 
the only study of the effect of a time-dependent ram pres¬ 
sure, as it would be the case for galaxies orbiting in clus¬ 
ters. Using a sticky particle code, they modelled the ram 
pressure as an additional variable force on the gas disc par¬ 
ticles. They observed a significant backfall of material that 
had been stripped during the time of maximum ram pres¬ 


sure but was not unbound . This is in agreement with the 
result of ISchulz fc Struckl (l200lll and iRoediger fc Hensleil 
(2005) that even for a constant ICM wind the unbinding 
of stripped mate rial t akes a few 100 Myr. The sticky par¬ 
ticle code of lVollmer et alJ 120010 cannot model hydrody¬ 
namical processes such as the continuous stripping found 
with hydro-codes because it cannot simulate hydrodynam¬ 
ical instabilities. Nonetheless, this model has been applied 
successfully to individual galaxies (e .g. IVollmer et al .11200(1 . 
I2001L Il 99 fll : IVollmerl l200(j l. iMarcolinjeteflTiKOO^I studied 
ram pressure stripping of discy dwarf galaxies in groups, 
where the ram pressure is lower than in clusters. They used 
a hydrodynamical 3D grid code and studied inclinations of 
0°, 45° and 90°. This group found that in dwarf galaxies 
the inclination angle only plays a role as long as the cen¬ 
tral disc pressure is comparable to the ram pressure. While 
for smaller ram pressures little gas is lost, for stronger ram 
pressures the complete disc is stripped. Even in the range 
where inclination plays a role, the 45° runs were similar 
to the fac e-on cases and only edge- on cases could retain 
more gas. IRoediger fc Henslerl (2QQ3) presented a compre¬ 
hensive parameter study of face-on stripping. They varied 
wind density and velocity and covered conditions from clus¬ 
ter centres to outskirts and galaxy groups. The ram pressure, 
Pram = PicmUicm: depends on, both, ICM density and veloc¬ 
ity, and usually it is assumed the stripping depends on p ram , 
and not on picm an d utcim individually. While t his was con¬ 
firmed in general bv lRoediger fc llerisleil i2005f) . they found 
that supersonic runs tend to suffer slightly less than cor¬ 
responding subsonic runs. This work used a 2D hydrody¬ 
namical grid code and was thus restricted to the face-on 
geometry. 

Even though ram pressure stripping has been studied 
by several groups, most simulations were performed for a 
face-on geometry, either because the code was 2D or be¬ 
cause this geometry is expected to show the strongest im¬ 
pact. Qualitative results for the influence of the inclination 
angle have been found, but they do not agree in all points. So 
here, we concentrate on the influence of the inclination an¬ 
gle in the case of massive galaxies and present 3D hydrody¬ 
namical simulations that were performed with the adaptive- 
mesh-refinement code FLASH. We study the amount of gas 
loss and discuss asymmetrical structures introduced in the 
remaining gas disc. We also present projected gas surface 
density maps which reveal interesting implications for the 
interpretation of observational data. 

In Sect. |5| we briefly introduce the code and describe 
our galaxy model and simulation parameters. In Sect. [3 
we discuss some analytical considerations, and in Sect.^lwe 
present our simulation results. Finally, in Sect.^lwe compare 
our results to previous work and discuss their implications. 


2 METHOD 

We study the motion of the galaxy through the cluster in 
the rest frame of the galaxy. Therefore, its motion translates 
into an ICM wind flowing past the galaxy. 


© 2005 ras, mnras ooo.mrHn 










































Ram pressure stripping of inclined disc galaxies 3 


2.1 Code 

The simulations we re performed with the FLASH code 
llFrvxell et alJ ? 2f)0(f i. a hydrodynamical adaptive-mesh- 
refinement code with a PPM hydro solver. All boundaries 
but the inflow boundary are open. The simulations pre¬ 
sented here are performed in 3D. We use a simulation 
box of size of (a; m m, £ ma x) x (f/min, |/max) X (Zmin, ^-max) — 
(—64.8 kpc, 64.8 kpc) 3 . For the low ram pressure runs a 

box Size of (iEmin, ^max) X (l/min, f/max) X (Zmin, £max) — 

(—97.2 kpc, 97.2 kpc) x (—64.8 kpc, 129.6 kpc) x 
(—97.2 kpc, 97.2 kpc) was used to prevent interference 
with the boundaries. The galactic centre is located at 
(a; g ai, J/gai, 2gai) = (0,0,0). The ICM wind is flowing along 
the p-axis into the positive direction. Most simulations 
were done using 4 refinement levels, resulting in an effective 
resolution of 500 pc and an effective number of grid cells of 
256 3 (384 3 for the low ram pressure runs). For comparison, 
we also run a few test simulations with 5 refinement levels 
(effective grid size of 512 3 , and effective resolution of 
250pc) (see Appendix^. The refinement criteria were the 
standard density and pressure gradient criteria. 

In the FLASH code, we used a mass scalar to “dye” 
the galactic gas. Thus, we can identify which fraction of 
gas inside a grid cell originated from the galactic disc (for 
a more detailed descriptio n of such “dyeing” techniques see 
iRoediger fc Hensieill2005lf . 


2.2 Model galaxy 

We model a massive spiral galaxy with a flat rotation curve 
at 200 kms -1 . It consists of a gas disc, a stellar disc, a stellar 
bulge and a dark matter (DM) halo. For the gravitational 
potential of the stellar disc, bulge and DM halo we use the 
following analytical descriptions: 

Stellar disc: Plummer-Ku zmin _ disc, se e 

iMivamoto fc Nagal (Il975h or iBinnev fc Tremain5 lll987fl . 
Such discs are characterised by their mass, M*, and radial 
and vertical scale lengths, a* and 6*, respectivel y. 

Stel lar bulge: Spherical Hernquist bulge (see iHernnnisl 
? 1 993h . In case of a spherical bulge, the gravitational poten¬ 
tial, <f>. depends on radius, r, as <f>(r) =- , Mbn ' gc , where 

x ' r_ r' r bulge 

Afbuige is the mass of the bulge, rb u i ge the scale radius and 
r the spherical radius. 

DM halo: The spherical model ofl Burkertl lll995h , including 
the self-scaling relations, i.e. the DM halo is characterised 
by the radial scale length rn M alone. For t he eq uati on of the 
analytical potential see also lMori fc Burkertl ( 200df) . 


The self-gravity of the gas is neglected as the gas contributes 
only a small fraction to the overall galactic mass. The galaxy 
model parameters are summarised in Table Q In order to 
prevent steep density gradients in the galactic plane and in 
the galactic centre, the gas disc is described by a softened 
exponential disc : 


P (R,Z) = 


Ma, 


27r a| as &gas 


0.5^ sech 


R 


sech 


^]. (!) 


The coordinates ( R,Z ) are galactic cylindrical coordinates. 
The radial and vertical scale lengths are a gas and 6 gas , 
respectively. For R > a gas and \Z\ > fe gaa , this density 
distribution converges towards the usual exponential disc 


Table 1. Galaxy model parameters. 



M* 

O 

1 

stellar disc 

a* 

4 kpc 


6* 

0.25 kpc 

bulge 

-Mbulge 

I 5 

o 

o 


^bulge 

0.4 kpc 

DM halo 

IDM 

23 kpc 


-Mgas 

I 5 

o 

o 

gas disc 

&gas 

7 kpc 


frgas 

0.4 kpc 


^rot 

200 kms -1 


p(R,Z) = A / ga ° —exp(—7?/a gas ) exp(-|Z|/b gas ). For the 

Z7ra gas°gas 

corresponding exponential disc, A/ gas is the total gas mass. 
We chose M gas such that in the outer regions the gas disc 
converges to an exponential gas disc with 10% of the stellar 
disc mass, i.e. M gas = 0.1 M*. Due to the reduction of the 
gas density in the central part of the softened exponential 
disc, the integrated gas mass amounts to 6 • 10 9 Mq. Given 
the density distribution in the disc, its pressure and temper¬ 
ature distribution are set such that hydrostatic equilibrium 
with the ICM is maintained in the direction perpendicu¬ 
lar to the disc plane. In radial direction, the disc’s rotation 
velocity is set so that the centrifugal force balances the grav¬ 
itational force and pressure gradients. We have cut the gas 
disc smoothly to a finite radius of 26 kpc to prevent interac¬ 
tion of stripped gas with the grid boundaries. This smooth 
cutting was achieved by multiplying the gas density distri¬ 
bution p(R,Z) with 0.5(1 + cos(n(R — 20 kpc)/6 kpc)) for 
20 kpc < R 26 kpc. Figure Q shows radial profiles of den¬ 
sity, surface density, pressure and rotation velocity for the 
initial model. 

We measure the inclination angle, i, of the galaxy as the 
angle between its rotation axis and the ICM wind direction. 
Thus 0° corresponds to to a face-on motion of the galaxy, 
and 90° to edge-on. 


2.3 ICM conditions 

In order to study the influence of the inclination angle, i, 
isolatedly, we use a constant ICM wind. The additional effect 
of a variable ICM wind, as one would expect for a cluster 
passage, will be studied in a forthcoming paper. 

The strength of the ICM wind is characterised by 
the ram pressure, p ra m = PicmUj CM , which depends 
on, both, ICM wind den sity, picm, and velocity, uicm- 
iRoediger fc Henslerl ( 2005 1 have shown that the ram pres¬ 
sure is the relevant parameter and not picM or uicm- There¬ 
fore, we use uicm = 800kms -1 for most simulations and 
vary the ram pressure by varying picm- We chose the ICM 
temperature such that the sound speed is 1000 kms' 1 . Thus, 
the standard vicm corresponds to a Mach number of 0.8. For 
comparison, a few runs with supersonic wicm = 2530kms -1 
have also been performed. The parameters for p ram , Picm, 
uicm and i are summarised in Table |3 We will refer to the 
three ram pressures used as low, medium and high. If not 
stated otherwise, we refer to runs with wicm = 800kms -1 . 

We start the simulation with the ICM at rest and then 


© 2005 ras, mnras ooo.rnrrn 

























4 E. Roediger and M. Briiggen 



0 10 20 30 

R / kpc 


Figure 1. Radial profiles of the density, p , pressure, p, and ro¬ 
tation velocity, u ro t, in the galactic plane for the initial model. 
Also the radial profile for the projected ISM surface density, E, 
is shown. 

Table 2. ICM wind parameters. Crosses indicate which ICM wind 
parameter - inclination combinations have been simulated. 


inclination i/° 

Pram/ PICm/ ^ICm/ 


(ergcm * * 3 * * * * * * * ) 

(gem 3 ) 

(km s 1 ) 

0 

30 

60 

75 

90 

6.4 • 10" 11 

10-26 

800 

X 

X 

X 


X 

6.4 • 10" 12 

to- 27 

800 

X 

X 

X 

X 

X 

<N 

1 

o 

CO 

O 

1 

to 

00 

2530 


X 


X 


6.4 • 10" 13 

O 

1 

to 

00 

800 

X 

X 

X 


X 


increase the inflow velocity over the first 50 Myr from zero 

to PICM- 

In order to include all stages of the ICM-ISM interac¬ 
tion, the overall runtime of our simulations is 1 Gyr. 

3 ANALYTICAL CONSIDERATIONS 

3.1 Different phases 

As described in previous works, the ICM-ISM interaction 

proceeds in multiple phases. The discussion of the phases in 

Sect. 0 made clear that the overall process is quite com¬ 

plex. Analytical estimates exist mainly for the instanta- 

neous stripping phase in face-on geometry (see Sect. m- 

iMarcolini et all 11200.31) estimated a stripping radius for the 


edge-on geometry by assuming that the instantaneous strip¬ 
ping can be neglected and the stripping radius will be set 
by the KH-instability. They estimated the radius up to 
which the galaxy’s gravity can suppress the KH-instability, 
and assumed that gas outside this radius will be stripped. 
They found that, both, the estimates for face-on and edge- 
on cases are rather similar. However, the timescale for the 
continuous stripping, which is the relevant process in the 
edge-on geometry, operat es on a much longer timescale. 
iRoediger fc Hensleil fcOOfill estimated that the mass loss rate 
during the continuous stripping phase is of the order of 
lM©/yr. 

In the following sections, we discuss the estimate for the 
stripping radius in the instantaneous stripping phase and try 
to derive an estimate for inclined geometries. 

3.2 Face-on geometry 

The usual way to estimate the amount of gas loss during 
the instantaneous stripping phase for galaxie s mov ing face- 
on thr ough the ICM follows the suggestion of lGunn fc Gotti 
lll972il . Here, one compares the gravitational restoring force 
per unit area and the ram pressure for each radius of the 
galaxy. At radii where the restoring force is larger, the gas 
can be retained, at radii where the ram pressure is larger, 
the gas will be stripped. The transition radius is called the 
stripping radius. In order to calculate the restoring force 
per unit area, it is assumed that the gas disc can be approx¬ 
imated by an infinitely thin disc. Then, the restoring force 
per unit area would be 

/ grav (R) = E gas (R) |§(R), (2) 

where >1? is the gravitational potential of the galaxy 
and S gas (-B) the gas surface density. As pointed out by 
e.g. iRoedieer fc Henslerl • 2IIII-V :. the gradient of dQ/dZ ex¬ 
actly in the disc plane is zero, so if d$/dZ exactly in the 
galactic plane is used, this estimate would not yield any 
restoring force. In order to obtain a reasonable estimate, 
one can use the maximum gradient ^|( R, Z ) behind the 
disc plane at a given R. In general, the resulting stripping 
radii and retained gas masses from numerical simulations 
agree well with the predictions from this simple estimate. 

3.3 Inclined cases 

Now we wish to move on from face-on geometries to inclined 
cases. As in face-on geometries, we assume that the gas disc 
can be described by an infinitely thin disc. For a given sur¬ 
face element dA at radius R in the galactic plane, the ram 
pressure force is picmUi CM cos id A, because we have to use 
the projected cross-section of this surface element. This force 
is axisymmetrical with respect to the galactic rotation axis 
for all i. We can also interpret this force in terms of an 
effective inclination-dependent ram pressure cos i Picm«icm- 
Clearly, the ram pressure force decreases with inclination. 
In addition to the reduced ram pressure force, also the cor¬ 
rect gravitational restoring force needs to be calculated. For 
a surface element dA in the galactic plane at galactic ra¬ 
dius R and azimuthal angle cj >, the restoring force is the 
projection of the sum of local gravitational force and cen¬ 
trifugal force onto wind direction. Since the centrifugal force 
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and the radial gravitational force in a disc galaxy balance, 
we can drop both terms from our analysis. Only the pro¬ 
jection of the gravitational force perpendicular to the disc 
plane onto wind direction remains. This is cos i • / gra v dA 
with /grav calculated as in Eqn. [5] Again, this force is ax- 
isymmetrical with respect to the galactic rotation axis. If 
now the ram pressure force and the restoring force are com¬ 
pared for all disc radii, the factor cosi drops out completely. 
Consequently, the stripping radius should be independent 
of inclination. One of the assumptions for this estimate is 
that the gas disc can be approximated by an infinitely thin 
disc. Clearly, this assumption is not valid for highly inclined 
cases, so the estimate is expected to fail for near-edge-on 
cases. However, for moderate inclinations it could give a use¬ 
ful approximation. For near-edge-on inclinations, we expect 
that the galaxy loses less gas. 


4 SIMULATION RESULTS 

The most striking result of our simulations is that with in¬ 
creasing inclination a stronger and stronger degree of asym¬ 
metry is introduced to the gas disc. This asymmetry makes 
simple concepts such as a stripping radius difficult to apply. 

4.1 Snapshots of the evolution 

To visualise the 3D data, we show slices through the simula¬ 
tion box and projected gas densities. Figures|2]and[]]display 
the local gas density in certain planes through the galaxy for 
an inclination of 30° (Fig. [SJ and 90° (Fig. [3>. Both plots 
represent runs with medium ram pressure. Snapshots at 3 
different times are shown. For the 30° case, the disc remains 
rather symmetrical, whereas especially early snapshots of 
the 90° case display a strong degree of asymmetry. 

In Figs. Q101 and 01 we show projected gas densities 
at several timesteps for the same runs as in Figs. 01 and 
0] but also for the corresponding 60° run. For all three 
runs, we show projections along the three coordinate axes. 
Only galactic gas (identified by its “colour”, see Sect. ITT! 
is shown. 

The snapshots reveal an interesting trend: At early 
times a higher inclination causes stronger asymmetry: one 
side (roughly the upstream side) is stripped more strongly 
than the other. Details of how the rotation of the disc af¬ 
fects the stripping are discussed below. However, even for the 
edge-on case, the asymmetry decreases with time and the re¬ 
maining gas disc becomes circular again. This is caused by 
the rotation of the gas disc. At a galactic radius of 10 kpc, 
the rotation time is roughly 300 Myr, so in the course of 
our simulation the disc rotates roughly 3 times. Therefore, 
for the case of a constant ICM wind, eventually each part 
of the disc passes the region of strongest stripping, and the 
disc becomes symmetrical again. 

4.2 Velocity information 

IVollmer et a,l] f 200 lh stressed that velocity information is 
crucial for the correct interpretation of HI observations of 
ram pressure stripped galaxies. Due to the fact that we 
have used a constant ICM wind, the velocity structure of 
the stripped material reflects mainly the superposition of 



0 10 20 30 

R / fcpc 


Figure 7. Evolution of radial profiles, for density, p, and rota¬ 
tion velocity, v TO t, in the galactic plane. Also the projected gas 
surface density, £, is shown. The profiles belong to the same run 
as in Fig. The thick line displays the mean profile averaged 
azimuthally over the disc, the thin lines display the minimum 
and maximum values for each radius. For further explanation see 
text, Sect. m The colours/linestyles code four different times, 
see legend. 


the galactic rotation and the acceleration by the ICM wind. 
The velocity structure will be more interesting for simula¬ 
tions with variable ICM winds. 


4.3 Radial profiles 

In Fig. |3 we display the evolution of radial profiles. To gen¬ 
erate this plot, we measured radial profiles along 12 radial 
lines that are separated by 30° in the galactic disc. The 
thick lines show the average of these profiles, the thin lines 
maximum and minimum value for each radius. The profiles 
reveal that the inner part of the gas disc remains symmet¬ 
rical, whereas the outer part becomes asymmetrical. 

The decrease of the central density is caused by the 
spatial discretisation of our grid. In grid-codes, advection 
of mass and momentum into non-axial directions is always 
accompanied by a numerical viscosity due to numerical dif¬ 
fusion. In the rotating gas disc, the gas has to move on cir¬ 
cular orbits in a Cartesian grid. This geometrical mismatch 
is strongest in the inner part of the galaxy. Below a cer¬ 
tain radius, the resolution of the rotational motion becomes 
insufficient and the numerical diffusion introduces radial ve¬ 
locities, which lead to a decrease of the central density. We 
checked in higher-resolution runs that this behaviour does 
not influence our results (see Appendix IXl. 
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Figure 2. Slices in y-^-plane at, x = 0. colour-coded density for different timesteps. For medium ram pressure, inclination of 30°. The 
ICM wind is flowing from left to right. Also the grid structure is shown: One square is the crossection of one block= 8 3 grid cells. These 
snapshots are taken from the run with 250 pc resolution. 



Figure 3. Slices in .r-y-plane at z = 0, colour-coded density for different timesteps. For medium ram pressure, inclination of 90°. The 
ICM wind is flowing from bottom to top. For more details see Fig.|3 


4.4 Mass loss / radius decrease 

4-4-1 Radius and mass measurement 

In the light of the previous plots, it is not clear how the 
radius of the gas disc should be defined. We decided to in¬ 
clude the asymmetry in our measurement. We scan along 
12 radial lines, separated by an angle of 30°, in the galactic 
plane. For each of these scan lines, we find the radius where 
the density drops below 10~ 26 gem -3 for the first time. This 
is the stripping radius for thi s scan lin e. T his density limit 
seems rather arbitrary, but lRoedigeJ (l2005h tested different 
options and found this one to be the most appropriate and 
representative one. We average the results of all scan lines 
to obtain a mean disc radius f? mea n, but we also compute 
the maximum and minimum radius f? max and Rmin- 


Two masses are of interest here: the mass Afdisc of 
gas inside a cylinder centred on the galaxy (radius=27kpc, 
thickness = ±5kpc), and the mass of gas bound to the 
galaxy, Mbnd- In order to exclude the ICM in Mdisc, we only 
consider gas with a temperature less than 10' K. The bound 
mass is given by the sum of the gas mass in all those grid 
cells where the total energy density etherm + eidn + e po t = 
p /(7 - 1 ) + 0.5 p(vl + vl + v 2 z ) + is negative. 

4-4-2 Radius and mass 

The evolution of disc mass and radius is shown in Fig. [H] The 
left, middle and right plot correspond to strong, medium 
and weak ram pressure, the inclinations are colour-coded. 
The reason why the initial mass appears less in the case 
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Figure 4. Projected gas densities along x -, z-, and 3 /-axis in the left, middle and right column, respectively. Each row is for the same 
timestep that is printed in the top right corners of each panel. For medium ram pressure and inclination of 30°. Contour levels are 
10 -5 ’ -4 '” 3 ’ -2 gem -2 . 


of the high ram pressure is that we have set the gas disc 
in pressure equilibrium with the ICM. For a higher picm 
but constant ICM temperature, the ICM pressure is higher, 
which causes the outer layers of the gas disc to have a higher 
pressure. As the density distribution is fixed, this means 
the outer layers have a higher temperature and, therefore, 
drop out of our calculation of Mdi SC du e to the temperature 
limit (see also discussion of this effect in lRoediger Rr. ffensletl 
l2005t i. In order to include the asymmetry of the disc, we plot 
Rmin(t) and i?max(t) as a function of time. In addition, we 
overplot the mean radius, i?mean(t). For strong and medium 
ram pressure, it is obvious that the degree of asymmetry in¬ 
creases with inclination. We have also plotted the quantities 


-Rmax Rm in and (Jimax Rulin')!Rmean in Fig.[l]tO quantify 
the degree of asymmetry. 

Several interesting statements can be made from the 
radius and mass plots: 

• Depending on ram pressure, we can strip the disc com¬ 
pletely, partially or only marginally. 

• The evolution of the retained mass and, particularly, the 
disc radius are remarkably similar for face-on, 30° and 60° 
runs. The runs for 0° and 30° are nearly indistinguishable, 
and in the 60° case only slightly less mass is lost. 

• For all but the edge-on runs, we can clearly disting uish 
the three phases discussed in lRoediger fc Hensler feoPS ^l and 
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Figure 6. Projected gas densities along x and y-axis. For medium ram pressure and inclination of 60°. See also Fig. PH 
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Figure 8. Radius and mass of the remaining gas disc for three ram pressures (left, middle and right plot correspond to strong, medium 
and weak ram pressure) and different inclinations (colour-coded). For the radius, we show the mean radius (thick solid line) and the 
minimum and maximum radius (thin lines of same colour). The area between the minimum and maximum radius curves is hatched for 
easier orientation. In addition to the colour, also the direction of hatching codes the inclination (e.g. horizontal hatching for edge-on). 
Thin dashed lines in mass plot are bound gas mass, thick solid lines are the mass in a fixed disc region (see text, Sect. n.'i. ii . Both, mass 
and radius do not drop immediately because we switch on the flow at the inflow boundary at t = 0, and then the flow needs a certain 
time to reach the galaxy. 





Figure 9. Demonstration of the asymmetry of the gas disc and its evolution. The asymmetry is shown by R m ax — /i m in in the top 
panels, and the same normalised to i? mea n in the bottom panels. For three ram pressures: left, middle and right plot correspond to strong, 
medium and weak ram pressure. In each plot, the results for different inclinations are shown (coded by colour/linestyle, see legend). 
Please note the different scales on the j/-axes of the bottom panels. 


in Sect.Q 

The first is the instantaneous stripping phase, during which 
the outer part of the gas disc is pushed in the downwind di¬ 
rection. This phase ends with the first minimum in -R mea n(t). 
The first minimum in Mdi S c(t) appears slightly later because 
the gas has to leave our defined “disc region” first before it 
is no longer considered as disc gas. 

The next phase is the intermediate phase, which is charac¬ 
terised by Mbnd(t) > Mdisc- 

Finally, the phase of continuous stripping follows. 

• The durations of, both, the instantaneous stripping 
phase and the intermediate phase are longer for weaker ram 


pressures. For a fixed p ram they are independent of inclina¬ 
tion. 

• For non-edge-on cases, the evolution of f? mea n after the 
instantaneous phase is slow. 

• For a given p ram , the mass loss rate in the continuous 
stripping phase seems to be independent of inclination. This 
feature is most relevant for medium ram pressures, where 
neither the mass loss is negligible (like for weak ram pres¬ 
sures) nor the disc is stripped completely (like for the strong 
ram pressure). 

The basic point in all these aspects is that the stripping 
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process is nearly independent of inclination as long as the 
inclination is not close to edge-on. 

The dependence of asymmetry on ram pressure and 
inclination is more complex. Moreover, the asymmetry de¬ 
creases with time. In general, the face-on cases show nearly 
no asymmetry. The absolute asymmetry (i.e. R m ax — 7?. m j n ) 
for the edge-on cases depends strongly on ram pressure, 
where higher ram pressures lead to significantly higher 
asymmetries. For the remaining inclinations, the absolute 
asymmetry is approximately independent of ram pressure. 
The dependence of the relative asymmetry (i.e. (R m s,x — 
Rmin)/Rmes.n) on ram pressure reflects the dependence of 
Rmean on p ram ■ The radial profiles in Fig. Q reveal that the 
asymmetry is restricted to the outer part of the gas disc. 
The inner part remains symmetrical. This is even true for 
the medium ram pressure edge-on case, as can be seen in 
Fig.EI Even for inclined cases, an inner part of the gas disc 
remains more or less untouched. 

The comparison between subsonic and supersonic runs 
with ide ntical ram pressure reveal s a similar result to the 
result of lRoediger fc Hensleri l2005lf . In the supersonic case, 
the galaxy can retain slightly more mass than in the sub¬ 
sonic case. The difference in retained masses increases with 
inclination. For highly inclined galaxies, also the asymmetry 
of the disc is temporarily larger in the supersonic case. 

4.5 Comparison analytical/numerical 

Due to the asymmetry of the disc and the ongoing evolution, 
it is not straightforward to compare the numerical results to 
the analytical predictions. A first question to answer is at 
which time we have to measure mass and radius in our sim¬ 
ulations to make a meaningful comparison to the analytical 
prediction. We take measurements at two timesteps: at the 
end of the instantaneous phase and at the end of the interme¬ 
diate phase. For the strong ram pressure, the intermediate 
phase is finished at t = 175 Myr, for the medium ram pres¬ 
sure at t — 400 Myr, and for the weak one at t = 800 Myr. 
Except for the edge-on cases, the stripping radius does not 
change during the intermediate phase. Therefore, the ra¬ 
dius measurements taken at the end of of the intermediate 
phase are also representative for the end of the instantaneous 
phase. The disc mass, however, changes during the interme¬ 
diate phase. As a further complication, R m ax(t) and Rmin{t) 
oscillate on rather short timescales. We average these quan¬ 
tities over a smoothing length of 100 Myr to derive typical 
values. 

In Fig. El we show the numerical results and the ana¬ 
lytical estimate for the stripping radius and the remaining 
mass in the disc. The analytical estimate is displayed by the 
dotted lines. The predicted stripping mass given in the lower 
panel of this plot is simply the disc mass inside the stripping 
radius. For the numerical stripping radius, we show R me an as 
diamonds and the range Rm in to 7? m ax as an error bar. For 
the mass remaining in the disc, we show the numerical re¬ 
sults at the end of the intermediate phase as diamonds. The 
disc mass at the end of the instantaneous phase is shown by 
stars. However, for the high ram pressure case this quantity 
cannot be derived. 

For the face-on cases (i = 0), the analytical estimate 
and the numerical result for the stripping radius agree well. 
The same is true for the disc mass at the end of the instanta- 



0 20 40 60 80 100 

i / degree 


Figure 10. The analytical estimate of stripping radius and strip¬ 
ping mass as a function of inclination for three different ram 
pressures (black=weak, red=medium, green=strong) is shown by 
dotted lines. Also the numerical results are shown. We measure 
the numerical stripping radius and mass at the end of the inter¬ 
mediate phase, see text for a more detailed discussion. For the 
numerical radius we give the mean radius as a diamond, and the 
range between the minimum and maximum measured radius as 
an error bar. See text for further explanation (Sect. 14.51 . The disc 
mass at the end of the intermediate phase is shown by diamonds. 
In addition, we plot the disc mass at the end of the instantaneous 
phase as stars. 

neous phase. The numerical values are just slightly smaller 
than the predicted ones. Here we confirm the result of previ¬ 
ous works. At the end of the intermediate phase the galaxy 
has already lost more than the gas outside -Rstrip- Therefore, 
the numerical values for the disc mass at the end of the in¬ 
termediate phase are significantly lower than the analytical 
prediction. 

For inclinations < 30°, neither stripping radius nor disc 
mass depend on inclination. For 30° < i < 60° the depen¬ 
dence on inclination is weak. Only for i > 60° the strip¬ 
ping radius and disc mass increase with inclination. There¬ 
fore, the simulation results confirm the prediction made in 
Sect. liOl 

The mass loss rate in the continuous stripping phase is 
remarkably similar in all cases and is about 1/ y r. This 
is similar to the results of lRoediger fc Hensleri ( 200? ). 

5 DISCUSSION 

5.1 The influence of the inclination 

According to our simulations, the inclination angle does not 
affect stripping significantly unless the galaxy moves close 
to edge-on. This result agrees with the tren ds found by 
\(£ ilis ct ah '2000i and iMarcolini et all 1200311 . In contrast 
to TouIikTetiidT ll200ll) . we do not find that only strict edge- 
on galaxies suffer less stripping, but we find that i affects 
the amount of mass loss for inclinations within ~ 30° of the 
edge-on geometry. Moreover, we have studied the stripping 
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process for longer runtimes and investigated a larger range 
of ram pressures. 

Our result disagrees with ISchulz fc Struck (200 11 who 
find that inclined galaxies are stripped on a longer timescale, 
also for inclinations around 45°. In contrast, in our simula¬ 
tions the phases of the ICM-ISM interaction are remark¬ 
ably independent of inclination. The d urati on of the p hases 
is com parable to what was found by iRoediaer fc Henslerl 

ft-U-1 

IVollmer et all (2001) have included a time-dependent 
ram pressure in their simulations, which seems to make the 
influence of inclination rather messy. They could not find 
a correlation between ram pressure, inclination angle and 
stripping radius (or remaining disc diameter), but they give 
an empirical fit for the remaining mass. We checked whether 
our simulation results for the disc masses follow their sug¬ 
gestion, but this is not the case. 

Based on their study of dwarf galaxies, iMarcolini et al, 
( 2003 1 argue that the inclination is only important when 
the ram pressure is comparable to the central pressure po 
of the galactic potential well. For p ra m > Po the galaxy will 
be stripped completely, independent of inclination, and for 
Pram < Po only very little gas is lost. For our galaxy, the 
central pressure is po « 10 -10 erg cm -3 . Our ram pressures 
are 6.4 • 10“ n ’~ 12 ’~ 13 erg cm -3 . If only the central galactic 
pressure is compared with p ram , a dependence on inclination 
should be prominent only in the case of our strongest ram 
pressure. In our simulations, however, at least the strong 
and the medium ram pressure case clearly depend on in¬ 
clination. We suggest that not only the central disc pres¬ 
sure is relevant for this comparison, but the pressure profile 
along the galactic plane. For our model galaxy, the pressure 
along the galactic plane drops to about 1CP 11 erg cm -3 at a 
radius of 5kpc and down to about 1CP 12 ergcm~ 3 at a ra¬ 
dius of 10 kpc. This explains why the inclination plays a role 
for the strong and the medium ram pressure. For the weak 
ram pressure, our galaxy loses little mass for all inclinations. 
However, even there, the influence of the inclination is the 
same. We conclude that the result of liVTarcolini et a.lJ (l2003ll 
cannot be applied to massive galaxies directly but needs to 
be improved in the sense that the inclination is important 
as long as the ram pressure is comparable to the pressure in 
the galactic plane. 


from the upwind edge has to go a longer way to reach a 
certain distance to the galaxy in p-direction, whereas the 
gas from the downwind edge has a “head start”. Moreover, 
in the direct vicinity of the galaxy, the ICM wind is not 
blowing straight in y-direction, but it is flowing smoothly 
past the galaxy. This introduces velocity components in zkx 
and ±z direction. Therefore, for some time during the initial 
phase the tail of stripped gas appears to be in the “wrong” 
direction. For such cases, not even radial velocity informa¬ 
tion would help to determine the true direction of motion 
because this behaviour is governed completely by velocities 
perpendicular to the line-of-sight. 

In the case of the map in Fig. El second line and second 
column, the reason is a different one. This galaxy is moving 
edge-on through the ICM, while we see it face-on. If we can 
determine the position of the galactic centre, e.g. from the 
stellar disc, we would see that the gas disc is stretched to¬ 
wards the top-right and compressed at the opposite side. But 
also this galaxy is moving towards the left. In this example, 
we can see the interaction between the galactic rotation and 
the ICM wind. In this map, the galaxy is rotating counter¬ 
clockwise. The side where the gas rotates parallel to the ICM 
wind is stripped more easily than the side where the gas ro¬ 
tates antiparallel to the wind. In later stages (further maps 
in the middle column of Fig. EJ, we can see that the tail of 
strippe d gas seems to be attached to the bottom side of the 
galaxy. IPhookun Kr Mundvl dl OQFill discussed the effects of 
combined stripping and rotation and concluded that in such 
a geometry the side where the gas rotates into the wind will 
be stripped more strongly because, as the relative velocity 
between the gas and the ICM is higher at that side, also 
the true ram pressure is higher at this side. Thus, the tail of 
stripped gas should be twisted towards this more strongly 
stripped side. This is just the opposite of what our simula¬ 
tions show. The relative velocity is higher at the side where 
the gas rotates into the wind, and also the ram pressure 
may be stronger there. But at this side the gas first has to 
be decelerated to zero velocity and then accelerated to get 
stripped. At the opposite side, the gas is already moving in 
wind direction and, hence, it seems easier to strip it there. 

In all cases discussed here, deep HI observations that 
reveal the low density extension of the tail or, for supersonic 
cases, the observation of a bow shock could help to determine 
the true direction of motion. 


5.2 What could be observed? 

A prediction for observable gas surface density maps is given 
in Figs. Hto© We want to emphasise one special aspect, 
namely the relation between the direction of the tail of 
stripped gas and the galaxy’s direction of motion. Consider, 
e.g., the following surface density maps (all at t = 200 Myr): 
Fig. El first line and first column; Fig.EJsecond line and sec¬ 
ond column; and Fig. El first colum. Also focus on contour 
lines of 10~ 4 gcm“ 2 and above, as it would be the case for 
“normal” HI observations. Now try to guess the galaxy’s di¬ 
rection of motion purely from these maps, from the direction 
of the tail. The answer is not that the galaxy is moving into 
the direction opposite to the tail, but in all cases the galaxy 
is moving towards the left. The reason for this strange ap¬ 
pearance is simple for the two cases where the disc is seen 
edge-on: Going back to the rest frame of the galaxy, the ICM 
wind is blowing the stripped gas towards the right. The gas 


5.3 Application to cluster galaxies 

The strongest limitation of our approach is the use of a 
constant ICM wind, which is not a good model for galaxies 
that traverse the cluster centre. Nonetheless, some of our 
results are applicable to real cluster galaxies. 

According to our simulations, the mass loss is similar 
for galaxies that do not move close to edge-on also in initial 
phase of stripping. Therefore, even if galaxies are exposed 
to strong ram pressures for a short period olny, only those 
galaxies that are close to edge-on can be protected against 
stripping by their inclination. Most galaxies should suffer 
severe stripping in cluster centres. 

Although the mass loss from galaxies in low ram 
pressure environments su ch as cluster outskirts is small, 
iRoediger fc llerisleil ifcOOfill found that still the gas discs are 
truncated at ~ 15 to 20 kpc. While this work was restricted 
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to face-on geometries, we found that this is also true for all 
non-close-edge-on geometries. Thus, ram pressure stripping 
of massive galaxies may also play a role in cluster outskirts. 
Observational evidence for the “pre-pr ocess in g” of galax¬ 
ies in low-density environments (see also lFuii tall2004 ) arises 
from the detection of HI deficient galaxies JSolanes et 33 
l200lll and passive spirals llGoto et a 3 l200.1l in cluster out¬ 
skirts. 

The duration of the intermediate phase, which is the 
phase where a substantial amount of gas that will be 
stripped eventually, is still bound to the galaxy, is inde¬ 
pendent of inclination. Therefore, we expect that if the 
ICM wind strength decreases during this phase, as it could 
be the case for galaxies passing cluster centres, a substan¬ 
tial amount of the stripped gas will still be bound to the 
ga laxy and will fall back to the disc. This effect was found 
bv IVollmer et alJ (l200lll in their sticky particle simulations. 
They also found that the amount of gas that falls back into 
the galaxy increases with inclination. It will be interesting 
to compare this feature with hydrodynamical simulations of 
a cluster passage. 

If the ICM wind decreases during the intermediate 
phase, it leaves behind a gas disc with a strong asymme¬ 
try in the outer parts. Differential rotation should smear 
out this asymmetry on a few galactic rotation timescales, 
resulting in a gas disc with a steepened gas density profile 
in the ou ter parts. A similar feat ure has been observed re¬ 
cently bv lKoonmann et alJ i2005h in the Ha discs of Virgo 
spirals. We suggest that this feature could be caused by the 
scenario we just described. 


6 SUMMARY 

We have presented 3D hydrodynamical simulations of ram 
pressure stripping of massive disc galaxies. We concentrated 
on the influence of the inclination angle between the galactic 
rotation axis and the galaxy’s direction of motion. 

We find that the inclination angle has a weak effect on 
the mass loss from the gas disc as long as the galaxy is not 
moving close to edge-on. We can explain this behaviour an¬ 
alytically by comparing the effective ram pressure with the 
component of the gravitational restoring f orce in wind di¬ 
rection. From simulations of dwarf galaxies. iMarcolini et alJ 
f 200fj l concluded that the inclination angle does not play 
a role as long as the ICM ram pressure is not comparable 
to the central disc pressure, po- For p ram < po hardly any 
gas would be lost, and for p ram > Po the gas disc would be 
stripped completely for all inclinations. We refine this result 
in the sense that for massive galaxies it is not sufficient to 
compare the central disc pressure with p ram , but the pres¬ 
sure profile along the galactic plane must be used for this 
comparison. This leads to a larger range of ram pressures 
for which the inclination angle plays a role. 

We have also studied the phases of the ICM-ISM in¬ 
teraction and found that for all non-edge-on geometries the 
stripping proceeds in a remarkably similar fashion: First, 
during the instantaneous stripping phase, the outer part of 
the gas disc is pushed towards the downstream direction. 
Then follows the intermediate phase which is needed to un¬ 
bind the stripped gas from the galactic potential. Finally, 


during the continuous stripping phase the gas disc contin¬ 
ues to lose gas at a small rate of ~ IMq yr” 1 . 

In addition to the mass loss, we have studied the asym¬ 
metrical structure introduced in the remaining gas discs of 
inclined galaxies. 

We have also presented projected gas density maps for 
several simulation runs. We have demonstrated that in mod¬ 
erately deep HI observations the tail of stripped gas does not 
necessarily point in a direction opposite to the galaxy’s di¬ 
rection of motion. Therefore, the observation of a galaxy’s 
gas tail may be misleading about the galaxy’s direction of 
motion. Deep HI observations or additional observations of, 
e.g., a bow shock may be essential. 
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APPENDIX A: RESOLUTION COMPARISON 

We have performed some of our simulation runs with two 
different resolutions to verify that our results are resolution 
independent. In Figs. IaTI and lA2l we repeat slices through 
the simulation box corresponding to Figs. [5] and [3] but for 
the lower resolution of 500 pc. In the higher resolution runs, 
the stripped gas fragments more than in the low resolution 
runs, but the overall structure is very similar for both reso¬ 
lutions. 

In Fig. |H we repeat Fig. Q but for the low resolution 
case. Here the numerical diffusion problem in the inner part 
as explained in Sect. l4.3l occurs in a larger region in the inner 
part. Aside from this, the characteristics of the profiles are 
very similar for the two resolutions. 

In Fig. we demonstrate that the disc radius and 
mass evolution are nearly independent of resolution. Even 
the mass loss rate in the continuous stripping phase is very 
similar. This is due to the fact that the mass loss is domi¬ 
nated by the largest KH-modes, which are also resolved in 
the lower resolution runs. 

We conclude that our simulations use a sufficient reso¬ 
lution. 
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Figure Al. Same as Fig. [21 but for lower resolution (Ax = 500 pc). 
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Figure A2. Same as FigE but for lower resolution run (Ax = 500pc). 
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Figure A3. Same as Fig. 0 but for low resolution run. 
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